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Summary 

High-explosive-driven  plasma  opening  switches  have 
been  modeled  in  one  dimension  using  the  Lagrangian  MHD 
code  RAVEN.  These  calculations  have  been  made  in  both 
cylindrical  and  planar  geometry.  Simple  compression 
can  account  for  observed  resistance  increases  at  early 
times  (time-of -flight  of  the  high-explosive  detonation 
products  across  the  plasma  conducting  channel).  Our 
results  suggest  that  some  improvements  in  switch  per¬ 
formance  might  be  achieved  through  a  judicious  choice 
of  gases  in  the  plasma  channel  and  by  lowering  the 
pressure  in  the  channel. 


Introduction 

To  optimize  the  operation  of  explosive-driven  mag¬ 
netic  flux  compression  generators  it  is  desirable  to 
have  a  low-inductance,  fast  opening  switch  that  can  in¬ 
terrupt  many  megamperes  of  current  through  a  ballast 
inductor  in  a  fraction  of  a  microsecond.  In  1977  Pav- 
lovskii  et  al.  reported  on  a  high  explosive  driven  op¬ 
ening  switch  with  which  they  achieved  a  resistance 
change  of  0.2  ohms  in  0.45  ps,  interrupting  a  current 
of  7.3  MA  and  switching  4  MA  into  a  30  nH  load. 

Experimental  efforts  to  examine  the  "Pavlovskii" 

switch  have  been  carried  out  by  Turman  and  Tucker  and 
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Turman,  Tucker  and  Skogmo  at  Sandia  National  Labora¬ 
tory  and  Goforth1*  at  the  Los  Alamos  National  Labora¬ 
tory.  When  working  at  current  densities  comparable  to 
or  higher  than  the  0.12  MA/cm  used  by  Pavlovskii  et 
al.  these  experimenters  have  seen  resistance  increases 
closer  to  50  mfl,  again  on  a  0.5  ps  time  scale. 

There  has  been  one  previous  attempt  to  simulate 
these  switches  computationally.  From  his  study  Baker3 
concludes  that  simple,  1-D  compression  of  the  current- 
carrying  plasma  channel,  the  mechanism  suggested  by 
Pavlovskii  et  al.  ,  cannot  account  for  the  observed 
resistance  increases. 


The  Model 

The  bulk  of  our  effort  has  also  been  a  one  dimen¬ 
sional  analysis  although  we  are  now  extending  our  ef¬ 
fort  to  two  dimensions.  Our  one-dimensional  calcula¬ 
tions  have  been  carried  out  using  the  1-D  MHD  code  RAV¬ 
EN6  in  both  cylindrical  and  planar  geometry.  RAVEN  is 
a  Lagrangian  code  that  uses  the  Braginskii '  formalisms 
for  electrical  and  thermal  conductivities. 

For  this  study  it  was  necessary  to  add  to  RAVEN  a 
high  explosive  equation  of  state.  The  internal  ener¬ 
gies  and  pressures  are  determined  from  a  burn  fraction 
which  is  determined  from  the  detonation  velocity.  For 
cyclotol,  which  we  used  in  calculating  cylindrical  ge¬ 
ometries,  we  used  an  internal  energy  of  9.2  x  10^6  ergs 
released  per  gram.  For  PBX  9501,  which  we  used  for 
planar  geometries  to  simulate  the  Goforth  experiments. 


we  used  an  internal  energy  of  1.02  x  10*1  ergs  released 
per  gram.  For  electrical  conductivity  of  the  HE  deton¬ 
ation  products  we  used  100  mho/m,  a  value  cited  by  Pav¬ 
lovskii  et  al.  For  thermal  conductivity  we  used  the 
Spitzer  formalism  assuming  the  detonation  products  to 
be  singly  ionized.  In  fact,  however,  varying  the  ther¬ 
mal  conductivity  by  a  factor  of  one  thousand  did  not 
significantly  effect  our  calculated  results. 

The  electrical  circuit  available  in  RAVEN  at  the 
time  of  this  study  is  shown  in  Figure  #1.  This  circuit 
will  not  permit  us  to  exactly  simulate  Pavlovskii *s  ex¬ 
periment.  We  can  introduce  an  external  load  on  the 
parallel  circuit  leg  (Lp),  but  this  inductance  will  not 

include  the  return  conductor  as  it  does  in  Pavlovskii *s 
experimental  set  up. 


TRANS  LINE 

GENERATOR  ♦ CABLES 


Fig.  1.  The  electrical  circuit  available  in 
RAVEN.  In  our  simulations  we  are  using  the 
generator  (L^  ),  the  ballast  inductor  (L^  = 
40  nH),  the  transmission  time  (L  *  30  nH),  the 
load  inductance  (Lj  =  1  nH),  and  the  parallel 
load  is  used  to  simulate  the  alternate  current 
path  so  that  Rp  is  dropped  to  zero  at  peak  com¬ 
pression. 


We  have  run  cases  in  both  cylindrical  and  planar 
geometry.  The  cylindrical  system  closely  approximates 
the  geometry  of  Pavlovskii.  As  is  shown  in  Figure  #2, 
the  cylinder  has  a  radius  of  12.5  cm.  The  plasma  col¬ 
umn  has  an  initial  radius  of  10  cm  and  may  expand  into 
a  0.5  cm  radius  gap.  The  central  conductor  has  a  radi¬ 
us  of  4  cm.  We  have  taken  the  cylinder  to  be  10  cm 
long.  Following  the  work  of  Baker3  we  note  that  after 
the  aluminum  foil  has  vaporized  the  aluminum  represents 
little  more  than  an  impurity  in  the  channel  which  we, 
therefore,  take  to  be  atmospheric  pressure  oxygen  at 
2.5  eV,  an  arbitrary  temperature  that  is  high  enough  so 
that  the  oxygen  will  conduct .  Our  planar  geometry 
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model  Is  taken  from  the  experimental  work  of  Goforth^.  than  that  measured  by  Goforth,  rather  surprising.  Our 
It  is  a  3cm  x  3mm  cavity  that  is  5  cm  in  length.  preliminary  analysis  indicated  that  these  switches 

would  be  stongly  limited  by  current  density  because  the 
higher  current,  and  subsequent  Joule  heating,  would  in¬ 
crease  the  ionization  level  of  the  plasma. 


Cell  Boundaries 


MHO  Tlm«  (pstc) 


'  central  conductor  -  plasma  channel 

- high  explosive  - external  insulator. 


Fig.  2.  Radii  of  each  zone  in  the  cylindrical 
calculation  as  function  of  time. 


Fig.  3.  Calculated  resistance  along  the  plasma 
channel  as  a  function  of  time  for  the  cylindri¬ 
cal  geometry. 


Results 

The  time  behavior  of  the  cylindrical  geometry  is 
shown  in  Figure  # 2.  In  these  calculations  we  hold  the 
outer  insulator  frozen  in  space.  As  a  consequence,  the 
results  of  the  calculations  after  the  initial  compres¬ 
sion  are  certainly  not  believable.  Figure  #3  shows  a 
resistance  increase  that  is  quite  comparable  to  the  in¬ 
itial  rises  reported  by  Turman  and  Tucker.  This  rise 
is  due  solely  to  the  compression  of  the  current  carry¬ 
ing  channel.  During  this  compression  the  temperature 
in  the  plasma  rises  sharply  to  nearly  14  eV  and  the 
pressure  reaches  just  over  1.8  megabars,  significantly 
higher  than  the  95-275  kilobar  pressures  in  the  HE 
zones . 

Our  resistance  drop  is  sharper  than  that  in  the  ex¬ 
perimental  results  of  Turman  and  Tucker.  Since  the  re¬ 
sistivity  of  the  ionized  plasma  is  inversely  propor¬ 
tional  to  the  temperature  to  the  3/2  power  this  sug¬ 
gests  that  our  calculated  temperature  Increase  may  oc¬ 
cur  too  rapidly.  This  would  not  be  surprising  since 
1-D  simulations  tend  to  overestimate  temperatures  at 
the  peak  of  compression  and  our  fixed  outer  insulator 
may  cause  an  overestimate  of  the  rate  of  compression. 

Our  calculated  secondary  resistance  increases  are 
smaller  than  the  first  and  this  does  not  agree  with  ex¬ 
periment.  In  addition,  these  secondary  increases  would 
probably  not  occur  if  the  outer  insulator  were  free  to 
move.  We  suspect  that  later  resistance  increases  are 
due  more  to  mixing  between  the  hot  plasma  and  relative¬ 
ly  cool  detonation  products  than  they  are  to  compres¬ 
sion.  To  verify  this  hypothesis,  we  are  presently 
carrying  out  a  series  of  2-D  calculations. 


Res i stance 


MHD  Time  (psec) 

Fig.  4.  Calculated  resistance  along  the  plasma 
channel  as  a  function  of  time  for  planar  geom¬ 
etry. 


In'  Figure  #4  we  show  the  results  for  the  planar  ge¬ 
ometry  with  a  current  density  of  0.22  MA/cm.  We  initi¬ 
ally  found  this  resistance  increase,  which  is  higher 


(3) 


Analysis  and  Conclusions 

If  we  are  dealing  with  an  ionized  plasma  (z,  >  1)  the 
resistivity,  q,  varies  as 


P  =  (y-l)U/V 

where  V  is  the  volume.  Simple  substitution  yields 


If  the  current  carrying  plasma  is  undergoing  an  adia¬ 
batic  compression  with  the  ideal  gas  y,  5/3,  then  the 

cross  sectional  area  will  also  vary  as  T  ^  2.  There¬ 
fore,  since  the  resistance,  R,  is  just 


2  Zeznz 

1+4 


3  (n.+n  )kT 
l  e 


(4) 


where  n^  is  the  total  ion  population  and  we  have  sub¬ 
stituted  in 


(1) 


P  =  (nA  +  ne)  kT 


(5) 


and,  the  length,  £,  does  not  change,  we  expect  the  re¬ 
sistance  to  be  constant  during  compression  once  the 
plasma  is  ionized. 

Figure  #5  shows  the  time  behavior  of  the  log  of  the 
pressure  as  a  function  of  density  for  the  planar  run  of 
Figure  #4.  The  line  plotted  during  the  initial  com¬ 
pression  shows  that  we  do  have  an  adiabatic  compression 
but  the  slope  is  not  1.67  but  closer  to  1.35.  This 
lower  value  of  y  indicates  some  sink  of  energy  not  in¬ 
cluded  in  the  above  analysis.  Additional  thought  indi¬ 
cates  that  this  sink  is  likely  to  be  ionization  of  the 
oxygen  plasma. 


Then,  in  terms  of  Z 


Y  =  1  + 


1  + 


2 


(l+Z)kT 


(6) 


In  the  limits  of  no  energy  in  ionization  (e  =0),  or 
very  high  temperature,  Eq.  (6)  limits  to  y  =  5/3  as  it 
should.  Between  these  limits  the  Saha  equation  must  be 
solved  based  on  the  density  and  temperature  in  the 
plasma.  We  have  solved  the  Saha  equation  iteratively 
for  values  of  constant  temperature  and  density  to  find 

values  of  n  and,  hence  Z,  in  oxygen  and  substituted 

these  values  into  Eq .  (6).  The  results  of  these  calcu¬ 
lations  are  shown  in  Figure  #6. 


Fig.  5.  Pressure  vs.  density.  For  a  Lagrangi- 
an  code  this  is  the  equivalent  of  pressure  vs 
(volume)”1.  Dashed  line  points  to  that  portion 
of  the  plot  that  was  traced  during  the  initial 
compression. 


Starting  Density=  1.00000e+14 


Armed  with  this  knowledge  we  can  estimate  y  based 
on  the  level  of  ionization.  The  internal  energy  can  be 
written 


Fig.  6.  Gamma  as  functions  of  temperature  in 
oxygen.  Each  line  is  a  factor  of  ten  in  dens¬ 
ity  higher  than  the  line  below  it.  The  bottom 
line  is  an  ion  density  of  10 14  cm  3,  the  top 
line  is  1023  cm  3. 


U  =  3/2P  +  l  ezqz 
z 


(2) 


where  ez  is  the  sum  of  the  ionization  energies  up  to 
level  z  and  r^,  is  the  number  of  ions  in  level  z.  We 
note,  however,  that  we  can  write  pressure  as 
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From  Figure  #6  we  reach  three  important  conclu¬ 
sions.  First,  it  is  possible  for  these  switches  to  be 

ionized  (Z  >  1)  and  still  have  y  <  1.67,  so  that  the 
resistance  can  increase,  even  if  weakly,  as  a  function 
of  temperature.  This  should  mean  that  these  switches 
can  be  used  at  higher  current  densities,  and  hold  these 
higher  current  densities  for  longer  times,  than  we  had 
initially  anticipated.  This  conclusion  appears  to  hold 
for  a  very  wide  range  of  temperatures.  This  tempera¬ 
ture  range  is,  however,  a  function  of  the  gas  in  the 
plasma.  Our  second  conclusion,  therefore,  is  that 
something  may  be  gained  by  carrying  out  an  analysis  for 
a  variety  of  plasma  gases  in  order  to  determine  which 
gas  will  provide  the  best  switch. 

Finally,  we  note  that  although  the  performance  of 
these  switches  in  the  compression  regime  should  be  rel¬ 
atively  independent  of  the  temperature  prior  to  com¬ 
pression,  over  a  wide  temperature  range,  Figure  #6  pre¬ 
dicts  a  strong  dependence  on  pressure  or  density.  This 
pressure  dependence  may  be  the  most  sensitive  parameter 
available  for  improving  the  performance  of  these 
switches . 
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